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VARIATIONAL PRINCIPLE AND STABILITY 
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Abstract. 

We prove a variational principle for stochastic Lagrangian Navier-Stokes 
trajectories on manifolds. We study the behaviour of such trajectories con- 
cerning stability as well as rotation between particles; the two-dimensional 
torus case is described in detail. 
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1. Introduction 

As discovered by V. I. Arnold ((2j) the motion of an incompressible non viscous 
fluid can be characterized as a geodesic on a group of diffeomorphisms. This point 
of view allows in particular to derive properties of the Lagrangian Euler flow, such 
as stability, through the study of the geometry of the group {[6]). 

When the fluid is viscous, namely for the Navier-Stokes equation, one can de- 
scribe the Lagrangian trajectories as realizations of a stochastic process and inter- 
pret the associated drift, solving Navier-Stokes, as an expectation over this process. 
This intrinsically probabilistic approach we follow here is inspired by [8] , [9] . Sim- 
ilar stochastic models are used for example in [4]. Then the trajectories remain, 
in an appropriate sense, geodesies and are almost sure solutions of a variational 
principle. This was shown in j3] for the two-dimensional torus. 

We prove a variational principle for the Lagrangian Navier-Stokes diffusions in 
a compact Riemannian manifold. Furthermore we study its stability properties. 
The behaviour of the trajectories depends on the intensity of the noise as well as 
on the metric of the underlying manifold. The example of the torus is studied in 
detail. Finally we describe the evolution in time of the rotation between stochastic 
Lagrangian particles. 
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Let (M, g) be a compact oriented Riemannian manifold without boundary. 
Recall that the Ito differential of an M-valued semimartingale Y is defined by 

(1.1) dYt = P(Y)^d( f P(Y);^odY, 



where 

(1.2) P{Y)^:Ty,M ^Ty,M 

is the parallel transport along t H> Y*. Alternatively, in local coordinates, 

il 
2 



(1.3) dYt = ( dY^ + -T]t,{Yt)dY^ ® dY,^ ] 



where r*^. are the Christoffel symbols of the Levi-Civita connection. 

If the semimartingale Yt has an absolutely continuous drift, we denote it by 
DYt dt: for every 1-form a E T{T*AI), the finite variation part of 

(1.4) / (aiYt), dYt) 

Jo 

is 

(1.5) / (aiYt), DYtdt) 

Jo 

Let G" , s > be the infinite dimensional group of homeomorphisms on M which 
belong to i/", the Sobolev space of order s. For s > ^ + 1: m = dimM, is 
a G°° Hilbert manifold. The volume preserving homeomorphism subgroup will be 
denoted by Gyi 

G^ = {<7eG^ : 5*M = m}, 
with fi the volume element associated to the Riemannian metric. We denote by ^^'^ 
(resp. ^#y) the Lie algebra of G" (resp. Gy). See [B] for example. 

On M we consider an incompressible Brownian flow Qu (t) G Gy with covariance 
a e r{TM TM) and time dependent drift u{t, •) e r{TM). We assume that for 
all X € M, a{x,x) = 2iyg^^{x) for some ly > 0. This means that 

(1.6) dgu{t){x) ® dgu{t){v) = a{gu{t){x),gu{t){y)) dt, 

(1.7) dgu{t){x) ® dgu{t){x) = 2j/g-i {gu{t){x)) dt, 

the drift oi gu{t){x) is absolutely continuous and satisfies Dgu{t){x) = u{t, gu(t)(x)). 
The generator of this process is 

d 

= i/A'' + — + 
ot 

where A'* is the horizontal Laplacian. The parameter v will be called the speed of 
the Brownian flow. 

If the time is indexed by [0,r] for some T > 0, we define the action functional 

by 



S{9u) = 



\Dgu{t){x)\\^ dx \ dt 



M 
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v{t,et{v)) 



2. The variational principle 

Define 

(2.1) Ji^ ^ {v e C\[0,T], ^^), v{0,-)=0, v{T,-)=0} 
Given v G J^f, consider the following ordinary differential equation 

detjv) 

(2.2) dt 

eo{v) = e 

where e is the identity of Gy. Since v is divergence free, e.{v) is a Gy -valued 
deterministic path. 

We denote by 3^ the set of continuous G^-valued semimartingales g{t) such that 
51(0) = e. Then for all v &Jif,we have et{v) o g^{t) G ^. 

Definition 2.1. Let J be a functional defined on and taking values in M. We 
define its left and right derivatives in the direction of /i(-) — e.(w), v £ at a 
process respectively, by 



(2.3) 



(i?L)hJ[5] = ^J[e.(e«)°5(-)]|.=o, 
iDR)hJ[g]^^J[9i-)oe.{ev)]U=o. 



A process g £ ^ wil be called a critical point of the functional J if 

(2.4) {DL)hJ[g] = {Dn)^,J[g], V/i - e{v), v e 

Theorem 2.2. Lef (i, a;) i^- u{t,x) be a smooth time- dependent divergence-free 
vector field on M, defined on [0,T] x M. Let gu(t) a stochastic Brownian flow with 
speed v > and drift u. The stochastic process gu{t) is a critical point of the energy 
functional S if and only if the vector field u{t) verifies the Navier-Stokes equation 

Ou 

(2.5) — + V„u = i^Du - Vp. 

For the construction of weak solutions of Navier-Stokes equations on Riemannian 
manifolds we refer to [?]■ 

Proof. Since the functional S is right invariant, it is enough to consider the left 
derivative. So we need to compute 

-I 



(2.6) 

We let 

(2.7) 

Then 

(2.8) 

which yields 
(2.9) /'(O) 



.oS{e.{ev){gu)). 
fie) ^ Sie.iev)igu)). 



/(e) = i^^(^E £ [\\De,{sv){gu)m^)\\ 



dt 



dx 




i{V,\e=oDet{ev) (g„(i)(a;)) , u(t,5„(t)(a;)))) dt 



dx. 
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We need to compute 

(2.10) Ve\e=oDet{ev){gu{t){x)). 
We have 

^ d , , ^ „ , det(sv) 
V.-|..oe.M=V.|..o^ 

= Vs\e=oev{t,et{sv)) 

= v(t,e). 

Together with ^^(O, •) =0, this imphes 

(2.11) ^\,=oet{ev){x)=v{t,x). 
Consequently 

(2.12) ^U=oet(ew) {9u{t){x)) = v {t, g^{t){x)) . 

By Ito equation, 

det{sv){gu{t){x)) 

(2.13) = {det{ev){-), dgu{t){x)) + ivdet(ef)(5„(t)(.T)) {dgu{t){x) ® dgu{t){x)) 
= {det{ev){-), dgu{t){x)) +v^et{ev){gu{t){x))dt. 

Here Aet(ew)(-) denotes the tension field of the map et{£v) : M — >■ M. This yields 
Detiev)igu{t){x)) = {det{ev)i-), Dgu{t){x)) + v^et{ev){gu{t){x)) 
+ ^v{t, et{ev){gu{t){x))) 
= {det{ev){-), u{t,gu{t){x)))+uAet{sv){gu{t){x)) 
+ ev{t,et{sv){gu{t){x))). 
Differentiating with respect to £ at £ = 0, we get 

Ve\e=oDet{ev){gu{t){x)) 

= {\/e\E=odet{ev){-), u{t, gu{t){x))) + e\e=o^et{ev){gu{t){x)) 

= ( V.-^|e=oet(£w)(-), w(i,S'„(t)(a;))) +z/n-^|e=oet(£w)(fl'«(t)(a;)) 

(2.15) \ ds / ds 

dv , , , , 
+ -^{t,gu{t){x)) 

dv 

= {W.v{t,-), u{t,guit)ix)))+unv{t,-){gu{t){x)) + —{t,gu{t){x)) 

dv 

= '^u{t,g^it){x))V{t, ■) + vOvit, ■){gu{t){x)) + —{t,gu{t){x)). 

We used the commutation formula V£|e=oA = D^, where □ = dd* + d*d is the 
damped Laplacian. Alternatively, 

(2.16) Dv = A^v + mc\v). 
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For a TAf-valued semimartingale J* which projects onto the M-valued semi- 
martingale Yt, we denote by ^Jt the Ito covariant derivative: 

(2.17) &Jt = P(Y)td {P{Y)T'Jt) ■ 
Then Ito equation yields 

(2.18) &uit,guit)ix)) ~ -^{t,g^{t){x))dt + S/ag^^t){x)U + i^A''u{t,g^{t){x))dt 
and 

(2.19) 9v{t,gu{t)ix)) ~ -^it,g^{t){x))dt + Vdg^(^t){.)V + i^A\{t,guit){x))dt. 
where the notation ~ means equal up to a martingale: 

P{gu{-));'&u{t,g^{t){x)) 

P{9u{-))T^ (^^{t,9u{t){x)) dt + ydg^(t)(x)U + z/A'^M(t, gu{t){x)) dt 

is a local martingale. 

On the other hand, denoting ut — u{t, gu{t){x)) and vt = v(t, gu(t){x)) we have 



{ut,vt)^ {Slut,vt)+ {ut,&vt)+ {^ut, &vt) 



(2.20) 



Let us denote by Dvt the drift of vt with respect to the damped connection on 
TM, whose geodesies are the Jacobi fields. It is known that, 



(2.21) 

and 

(2.22) 



(^Dut ~ iyRic^{ut)j dt is the drift of &ut 
(^Dvt - iyRic\vt)j dt is the drift of &vt. 



As can be seen from (|2.15p . (I2.19P and (I2.22p . the drift Dvt commutes with the 
derivative with respect to a parameter, so it satisfies 



(2.23) 



Dvt = Ve\e=oDet{ev){gu{t){x)). 



Taking the expectation in ((^^ and using (g^Sl), and (1^^ . we get by 

removing the martingale parts 



E[{ut,vt)] = E 



(2.24) +E 



E 



{-g^{t,gu{t){x))+Vu,u + i^A''u{t,gu{t){x)), vt) dt 



{ut, ye\e=oDetiev){gu{t){x)) - i^Ric\vt)) dt 



2v I tT{V.u, V.v) {t,gu{t){x))dt 



Then using the facts that vt = 0, together with 
(2.25) {ut,RicHvt)) = {RicHut),vt) 
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and p.l6p . we get 



E 



(wt, Ve\e=oDet{ev){gu{t){x))) dt 
^ ,du 



dt 



(2.26) 



Integrating with respect to x yields 
(2.27) 
/'(O) 



(t,gu{t){x)) + VufU + vUu{t,gu{t){x)), Vt) dt 



2v \ tr(V.Ut, V .vt)(t,g^(t)(x))dt 



-E 



E 



\Jm 

T 



+ V„ + i^n U (t,.g„(t)(x)), v{t,gu{t){x))) dx \ dt 



2v 



dt 



tr(V.u, \I .v) {t,gu{t){x))dx] dt 



10 \J M 

Now we use the fact that gu{t){-) is volume preserving 
/'(O) 

- rT , r / / / a 



(2.28) 



-E 



■E 



\Jm 

T 



+ V„ + iyD]u] {t, x), v{t, x)) dx] dt 



2v 



dt 



tr(V.u, V.-y) {t,x)dx dt 



/o \^A/ 

Since M is compact and orientable, an integration by parts gives 



(2.29) 



/ tr(V.u, .v) {t,x) dx ~ — / {Du,v) {t, x) dx. 
Jm Jm 



Replacing in (|2.28p we get 



(2.30) /'(O) = -E 



V"' M 



-Q^ + — lyD ) It ) {t,x), v{t, x) I dx ] dt 



The process gu{t) is a critical point of the energy functional S if and only if /'(O) 
0, which by equation (|2.30p is equivalent to 



(2.31) 



d 

— + V„ - ) It = -Wp 



for some function p on [0,T] x Af. This achieves the proof. □ 

3. A MARTINGALE CHARACTERIZATION FOR SOLUTIONS OF NAVIER-StOKES 

EQUATIONS 

In this section, to simplify the equations, we assume the pressure to be constant. 
The pressure will not be present, in any case, in the weak version of the formulae 
we derive. 
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We seek to obtain a formula for the drift of the covariant derivative with respect 
to a parameter of a family of Navier-Stokes solutions, extending the well-known 
Jacobi equation. 

Consider a family of diffusions g", a G R, satisfying 

(3.1) g"(0)=^(a) 

where ip -.R ^ AI is a smooth path on M, and solution to the Ito SDE 

(3.2) dg''{t) = u{t,g"{t))dt + a{g^{t))dBt 
where u solves 

(3.3) dtu + V„M + vDu = 0, 

Bt = {Bf)e>o is a family of real Brownian motions, a — {ae)e>o, and for all £ > 0, 
a I is a vector field on M . We furthermore assume that 

(3.4) aa* = vg-^ 

where g is the Riemannian metric on M . 

We denote by = Og'^it) = u{t,g°'{t)) the drift of g". We denote by ^"Jt 
the vertical part of the Ito differential (with respect to V^) of a TM- valued semi- 
martingale Jt- It is known that 

(3.5) ^Vt = &Jt + ^RiJt,dXt)dXt 

where Xt — 7r( Jj) and R is the curvature tensor. If Jt has an absolutely continuous 
drift D'^Jt, then the finite variation part of ^'^Jt is D'^ Jt dt. 
From the Ito equation 

(3.6) ~ dtu{t,g"{t))dt + \7dg'=-{t)U + iynu{t,g°'{t))dt 
we deduce that the drift of S''^uf is 

(3.7) = dtuit,g''{t)) + Vufu + lyDuit, g''{t)) = 0. 
^From [T] Theorem 4.5, we have formally 

(3.8) 

=Vo.^uf + i?(d(7«(t), a„5"(<)X 

+ R{dg^{t),d^g''{t))i^u^ - ,yd*R{d^g"{t))u^ + ^R{dg"{t),9d^g"{t))u?. 

Using (13.51) . we obtain 
(3.9) 

=Va^Sr + i?(d5"(i),aag"(t))u^-i.Va„g.(t) Ric««)dt 

+ R{dg"{t), do.g"{t))&u? - iyd*R{d^g"{t))u? dt + ]^R{dg^{t), &dc.g"{t))ut . 

Removing the martingale part we obtain the drift 
(3.10) 

D-Vc^uf =V^D-uf + R{uf,dc,g^{t))u? - vVe^g^(t) Ric*(wn 

+ 2i.tr i?(-,a,g"(i))V.u," - vd*R{do,g''{t))u'^ + J.tri?(a(-), Va„g=(t)a(-))Mt". 
Now since Duf = we finally get 
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Proposition 3.1. The drift of the covariant derivative with respect to a of the 
family {uf)aeR of Navier-Stokes solutions is given by 

(3.11) 

D'Vo^u? =R{u?,dc.g"{t))u? - y^a^g^it) Ric«(?/?) 

+ 2i. tri?(., 9„5"(i))V.M," - vd*R{do.g''{t))ut + J^tri?(a(-), Va„g.(,)a(.)X. 

This formula extends the weh known corresponding (Jacobi) equation for the 
variation of geodesies. 

4. The two-dimensional torus endowed with the Euclidean distance 

We study the evolution in time of the distance between two particles in the 
two dimensional torus. Notice that, in order to interpret the diffusion processes as 
a solution of the variational principle described in section 2, there is no canonical 
choice for the Brownian motion, as far as it corresponds to the same generator. We 
make here a particular choice. 

On the two-dimensional torus T — M/27rZ x R/27rZ we consider the following 
vector fields 

Ak{e) ^ (fc2,-fci)cosfc.0, Bk{e) = {k2,-ki)smk.0 
and the Brownian motion 

(4.1) dWit) = ^kV^Wxk + Bkdyk) 

where Xk,yk are independent copies of real Brownian motions. We assume that 
Sfc I^P-^fe ^ necessary and sufficient condition for the Brownian flow to be 

defined in L^(T). Furthermore we consider Xk — A(|fc|) to be nonzero for a equal 
number of ki and fc2 components. In this case the generator of the process is equal 
to 

d 
ot 

with 2C = -^fc (c.f.[3] Theorem 2.2). We shall assume C to be equal to one. 
Let us take two Lagrangian stochastic trajectories starting from different diffeo- 
morphisms and let us write 

(4.2) dgt = iodW{t)) + u{t, gt)dt, dgt = (odW{t)) -t- u{t, gt)dt 
with 

.90 = 0, 9q = iI^, (I'T^i^ 
We consider the distance of the particles defined by 

p2(0,v)= / \m-m\'de. 

where d0 stands for the normalized Lebesgue measure on the torus. 
We let pt = p{9t,gt) and T(g,g) = inf{t > : pt = 0}. 

Lemma 4.1. The stopping time T{g,g) is infinite. 
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Proof. By uniqueness of the solution of the sde for gt we can let for all t > 
9t{(^) = 9t{{4'~^ ° ^){()))- Since gt, v' and V are diffeomorphisms, if ip{6) ^ 'ip{6) 
then gt{0)^gt{{r'oij){e)). 

Since </> 7^ ^, the set {9 € T, gt{d) ^ 9t{S)} has positive measure and this implies 
that pt > 0, which in turn implies that T{g, g) is infinite. □ 

Denote by Lt{6) the local time of the process \gt{0) — gt{6)\ when {gt{6) , gt{9)) 
reaches the cutlocus of T. By Ito calculus we have 



dpt =—y^^^kVy {gt - 9t,{Ak{gt) - Migt)) dxk{t) + [Bkigt) - Bk{gt))dyk{t))j 

1 If 

H {9t-gt,u{t,gt)-u{t,gt))T dt / \gt - gt\iO)dLt{9) 

Pt Pt Jt 

+ ^ E^fe^ (\\M9t) - Ak{gt)\\l + WMgt) - Bk{gt)\\l) dt 

~ Y^^^l^ iS^t - gt, Ak{gt) - Ak{gt))l + {gt - gt, Bk{gt) - Bk{gt))l^ dt 



where (•, •)t and || • ||t denote, resp., the L inner product and norm. We let 

(4.3) Su{t) = - {u{t, gt) - u{t, gt)) . 

Pt 

We have 

(4.4) Mgt) - Mgt) = -2sin ^:%±M sin ( ti^) k\ 



iAK\ T3 ( \ T}(~\ o k-{gt+gt) . (k-{gt-gt)\,_i_ 
(4.5) Bk{gt) - Bk{gt) = 2cos — ^— sm I — ^— I k-^, 

where we have noted fc-"- = {k2, —ki). Then, for fc 7^ we let 



k 1 
(4.6) rife = — , and ng{t) = —{gt- gt). 

\k\ Pt 



This yields 



,2 ■ k-{gt+gt) sin (k-{gt-gt) 



(4.7) Akigt) - Ak{gt) = -2|fc|^p,sin |fc|p« \ 2 ' ""^^^ 



(4.8) BM - B,{gt) = 2|fc| V.(cos ( '-^^ ) 
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With these notations we get 
dpt 



sin [ k ■ {gt{e) - UO)) 



\k\pt 



= Pt^^Afc|fc|2^2(nfc^ •ng(t,0)) 

X (- sin'^ • + d.,it) + cos ^ • dy^jt) ) dO 

+ Pt{ng{t),Su{t))jdt-pt [ \ng{t,e)\-dLt{e) 

Jt Pt 



sin fk - {gt - gt) 



\k\pt 



dt 



dO dt 



+ 2vptY,>i\k\ 

k 

-2vptY,W 

And finahy: 

Proposition 4.2. T/ie /id equation for the distance pt between the diffeomorphisms 
gt and gt is given by 

(4.9) dpt = Pt {(Ttdzt + btdt+ {ng{t), 6u{t))j. dt - dat) 

where Zt is a real valued Brownian motion, at > is given by 
(4.10) 



de dt. 



X l^J (nfcj. • ng(t,^))sin 

k 

the process bt satisfies 



■{9t{0)+gtm \ sin f k-{gt{0)-~gt{0)) \ 
2 J \k\pt 



k ■ {gt{0) + ~gt{e))\ ^ f k ■ {gt{9) - ^9)) 
\k\pt 



de 



de 



(4.11) bt + -al=2vptY,>l\k\' 



sin (k - {gt - gt) 



\k\pt 



dt 



and at is defined by 

(4.12) ao = 0, dat= [ \ng{t,e)\—dLt{e). 

Jt Pt 
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Pt = Pto exp ( / asdzs+ I [bs~ -a^ + {ng{s), 6u{s)),^ ] ds - {at - ajj ) . 



Pt{ng ■ rik) 



So we have for all < to < t, 
(4.13) 

■ 1 
2 

Let 
(4.14) 

Notice that 

Lemma 4.3. We have 

1r^ ^2 ^4„\^\2|,|4 f ,2 sin^ ( k ■ {gt{e) - ~gt{e)) \ 

(4.15) at<Av^X,\k\ ^ J 

and 

(4.16) . > 2.Ea^i*i'X(«« «.)^^« (ttoMl) 



"* ^ / • I * I r7 I — I i±i n \ \ \ 

d0, 

k 

in particular 6* > 0. 

Let R> 0. Assuming that Afc = for all k such that \k\ > R then on 



{a;|V^eT, \g,{e) - ~gt{e)\ < ^} 



we have 

bt - > 0. 
Proof. Using Cauchy Schwartz inequality, 

a, <4^A,|fc| uj^\{ng-n,.)\sm ^ j ^ ^ j de 

' A '^"^ ■ "'^"^'^ I 2 ) 

+ 4i:A^.|.^/|(n, .n..)|cos^ (M,*(^I±m) ^ (M^^IzM))) 

x/L-n..)|^f^i^^(^^V^ 

' ■ \ /„p ■ ~ " " ifcim \ y. I 

k 



\k\pt 

=4.1:^11^-1' (Xic......)|igi( ^-'^'W-^'"» ; 

V VT Pt \k\Pt V 2 ; , 



12 



M. ARNAUDON AND A. B. CRUZEIRO 



On the other hand, 

_ o,,y^ \2i,,i4 /■ sin^ f k ■ {gt{9) ~ gt{0)) 

k 

SO that using the bound 



for (7^ yields 



where we used 



Since depends only on |fc|, we have \k — Aj,-l for all k. Then putting together 
the terms corresponding to k and k^ we obtain 



h + U^uY^\i\k\' 



2 

k 



sin^ ( k ■ {gt{9) ~ U0)) \ ^ ^ ^ fc^ ■ (.9,(0) - g,(g)) ^ 
and this yields using the bound for as well as 5^ + (5^^ = 1 

b,-y,>uY.xi\k\^ 

k 

.2^/^sin2 ( h-{g^{9)-~gm \ sin^ ^fc-^ . (^^(g) _ ^^^^^^ 



d9 



k 

V /^A^ A2 ^ ^in^ / |fc||.g,~.g,|(g) ^ sin^ / \k\\gt ^ ~gtm \ 
Assuming that = whenever |fc| > i? then on 



d9 



{c.|V0eT, \gt{9) - ~gm < ^} 



the functions inside the integral are nonegative, consequently 

1 



bt - > 0. 



Define 

(4.17) £{x) = ^ for x^O, £{0) = 1. 

X 

From Lemma 14.31 we easily get the following result. 



□ 
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Proposition 4.4. Let R > 1. Then on 

7^^/2l 



uj\y0eT, \gt{e) ~ ~gtie)\ < 



R 



letting 

" |4 



« = (^) E ^im 

^ ^ fc <fl 



we have, 

(4.18) dpt > pt (atdzt - ||(5u(<)||Trfi- / 6*) 1—^X4(6*) + 

V Jt Pt 

Moreover assuming that Xk — whenever \k\ > R, then letting 
o V —1 — 



|f| = - 

|fe|<_R 



on 



2R. 

(4.19) dpt > Pt (^cTtdzt + ^cr^dt- \\6u{t)\\T dt + c'j^ dt 

Proof. If 134(6*) - gti9)\ < then for all k such that \k\ < R, 

R 

k-{gt{e)~~gtm 



and this implies 



sin^ ^fc.(gt(0)-gt(0))^ ^ 1 .2 ^ ^ 



— ] {nu ■ ng) 



\k?pl V 2 ; ^ 4^ VV2 

So with KWf we get 



4 VV2 



2 



fc|<i?. 

k\<R. 



(again we used Jj{ng ■ n^)^ d0 + /^.(ng • n^.^)^ d^? = 1.). This establishes (|4.18p . 
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Next if \gtiO) — 5t(^)| < ttt; then from the calculation in the proof of Lemma l473l 



fe|<i?, 

^ / x2 \f s^"' fx M9t~~9tm \ sin' / \k\\9t-~9tm \ , 



IfcKi? 



M ^ D "^TT 



this estabhshes (|4.19p . □ 
Theorem 4.5. Lei t > 0. i? > 1 anrf 

// we assume the initial conditions for the distance and the l? norm of the initial 
velocity related as c = po — 2||uo||t > 0, and suppose that u = 0, then on fit, 

(4.20) ys<t, > e-''o'"=''^=+=«* ^po - 2||uo||t^ e-^°'"'-''^'-('=«+5)"ds 

as long as the right hand side stays positive. 

On the other hand if we assume that there exist constants ci , C2 > such that 
for alio &T and s € [0,i] , 

(4.21) |Vw(i,6l)| < cie"^^*, 
then on ^lt, Vs < 

(4.22) ps > Po exp ( / as dzg + c^t - — (l - e" 

\Jo C2 

Proof. Assume that po — 2||7io||T > 0. From inequality (|4.19p we have on fit, for 
s<t, 

(4.23) dps >Ps (^(Jsdzs + {c'j^ + ^a^,) ds^ ~ 2\\u{s, ds. 

Using the fact that u{t, .) satisfies Navier-Stokes equation together with Poincare 
inequality, 

^\\uis,.)\\l^^2i,\\\7u{s,.)\\l 
ds 



<-'^Ms,.)\\l 



Therefore we have 



\\uis,.)\\T<e-^'\\uo\ 



We obtain 

(4.24) dps > Ps (^cTsdzs + {c'ji + ^crf)ds^ - 2e""^^||wo||Trfs. 

From this comparison theorem for solution of sde's yields (|4.20p . 
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Now assume (|4.2ip . To prove (|4.22p we start with ()4.19p . and remark that 

Pw(OI|t < sup |Vu(i, 6')|. Then with the bound on \/u{t,9) we have 
def 

dpt > Pt ^cTtdzt + -CTj dt — cie~'^^* dt + Cj^ dt^ . 
Integrating the right hand side between to and t gives the resuh. □ 

Remark 4.6. The bound (|4.2ip is satisfied for instance for solutions u{t, •) of the 
form e-'^l'^I'Mfc. 

Also notice that, by the expression of the constant c'j^, the stochastic Lagrangian 
trajectories for a fluid with a given viscosity constant tend to get apart faster when 
the higher Fourier modes (and therefore the smaller lenght scales) are randomly 
excited. 



5. The two-dimensional torus endowed with the extrinsic distance 

It seems difhcult to deal with the local time term of Proposition To circum- 
vent this problem we propose to endow the torus T with a distance pr equivalent 
to the one of section |4j but such that p^ is smooth on T x T. Then we will see that 
when the assumptions of Theorem 14.51 are not fulfilled, then the behaviour of the 
distance of two diffeomorphisms can be completely different even if their distance 
is small. So the uniform control of the distance in Theorem l4.5l looks as a necessary 
condition for an exponential growth of the distance. 

The map 

ttZ X R/27rZ [0, 2] 
(6*1,612) ^ 2 



defines a distance on the circle R/2ttZ: it is the extrinsic distance on the circle 
embedded in the plane. From this distance we can define the product distance on 
the torus T. 



2 I '^1 111 „;„2 



1/2 



Pt((6'i, ^2), {0'i,0'2)) = 2 I sm" I I + sm- 

Note 

pU{0i,d2), {d[,e!,)) - 2 (2 - cos{9[ - 9,) - cos{9', ~ ^2)) ■ 

The distance p| is smooth on T x T. Now let cj) and be two diffeomorphisms on 
the torus T. We define the distance p{(j), ■0) with the formula 

pH(b,^)= f pum,m)d9 

= 2 J (2 - cos(0i(6') - ^1(6')) - cos(02(6i) - 0^(61))) d9 

=4] (sin^ (^m_m'^ + ,i,2 i^m_m^^i 

Now let 

Pt = p{9t,gt)- 
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^From the smoothness of pj, the formula for pt does not involve a local time. More 
precisely, letting 

Sg = gtie)-gtie), 
Scosk ■ g — cos A: • gt{(f) — cosfc • 54(6*), 
b sin k ■ g — sin k ■ gt{9) ~ sin k ■ gt{9), 
s\a.5g= {s\n{6gt)i{0), s\n{5gt)2{0)), 
Su = {u{t,gt) - u{t,gt)) 

we get from Ito calculus 

\ /sin(55 fScosk-g Ssink-g 

dpt = Pt; Afe ( ,(k2,-ki)\ dxk^ dyk 

^ \ Pt \ Pt Pt 

I sin 6n Su\ 
+ Pt( — ) dt 

\ Pt Pt 1 1 

This clearly has the form 

dpt = Pt {crt dzt + bt dt) 

where a, and bt are bounded processes and Zt is a real- valued Brownian motion. 
However it can happen that the drift is negative even if pt is small, as the following 
example shows. 

Example 5.1. Let a > small and e > satisfying e « a. Take (p = id and 
assume that there exist two subsets Ei and E2 of T such that Ei C E2, Ei has 
measure a, E2 has measure a + e, tp^O) = for all 9 G T\E2 and il){0) = {Oi+n, 62) 
for all 6 G El. Since s can be as small as we want, we have 

/Oo-4a, {smdg)o^O, {dgo)2^0, {6smk ■ g)o '=^0, 

on T\iJ2, {6 cosk ■ g)o = 0, 
on El, {S cosk ■ g)o = —2 if Aiiis odd, ((5 cos fc • 5)0 = if fciis even, 
so at time t = 0, 



dpt - -f I E ^^^i I 

\ki odd 



To construct a diffeomorphism like ip, one can cut an annulus Ei of width — in 

T and rotate it by tt. This yields a one to one map on T. Then smoothen it around 

the boundary of the annulus to get ip. The set E2 can be taken as an annulus of 

, , a + e ^ 
width — - — containing Ei. 

ZTT 
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6. Distance and rotation processes of two particles on a general 

riemannian manifold 

6.1. Distance of two particles. 

Let Bt = {Bf)e>o be a family of independent real Brownian motions, a = 
{o'e)e>o, with, for all i>0,aea divergence free vector field on M. We furthermore 
assume that 

(6.1) a{x)a*{y) = a{x,y). 
In particular 

(6.2) aix)a*{x) =2i^g-^{x). 
We lot (fjipG Gy. In this section wc assume that 

(6.3) dgt{x) = a{gt{x)) dBt + u{t, gt{x)) dt, go = (p 
and 

(6.4) dgt{x) = a{gt{x)) dBt + u{t, gt{x)) dt, go = tp 
For simplicity we let Xt = gt{x), yt = 9t{x) and 

Pt{x) = PM{xt,yt) 

For x,y £ M such that y does not belong to the cutlocus of x, we let a Ja{x, y) 
be the minimal geodesic in time 1 from x to y (70(2;,?/) — x, ji{x,y) = y)). For 
a e [0, 1] we let Ja = T^a the tangent map to 70. In other words, for v G T^M and 
w G TyM, Jaiv^w) is the value at time a of the Jacobi field along 7. which takes 
the values v at time and w at time 1. 

We first consider the case where yt does not belong to the cutlocus of Xt- We 
note To = Ta{t) = 7a (a;*, yt) and ja{t) = Ja{xt, yt)- 

Letting P{'ya)t be the parallel transport along Ja{t), we have for the Ito covariant 
differential 

^7a(t) := P{la)td{P{7a)t'ia{t)) 

On the other hand the Ito differential dja {t) satisfies 

djait) = Ja{dxt,dyt) + ^ {'V(^dxt,dyt)Ja-) {dxt,dyt). 

So we get 

(6.5) ^ja{t) = Vj„(da.,,dj/,)7a + Vi(v^^^^^^^^^^_^^(^^^ ,^^^)7a + iVd^„(t)-Vd^„(t)7a(i)- 

Let e{t) e Tx^M be the unit vector satisfying To{t) — pt{x)e(t). For £ > wc let 
a !->■ J^{t,x) be the Jacobi field such that jQ{t,x) = ae{gt{x)), Jf(i) = ae{gt{x)). 
Moreover we assume that Vj«(j 3.) Jo(t, a;) = and ^ jt(t,x)Ji{'^'X) = 0. 
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With these notations, equation (|6.5|) rewrites as 



l>0 ° £>0 



ja{dxt,dyt) + ^y^yjiW jeTgdt. 



2 

£>0 



We have 



.1 ^ 1/2N 

t^Pt(a;) I ( / {Ta{t),Ta(t)) da 



1 



2pt{x) 
1 



3d(||rof).d(||Toir) 



8/9t(a:) 

^ U^{t,x),et{x)\ dBf + (jo{u{t,gt{x)),u{t,gt{x))),et{x) 



Note 



/ {V jiW j/Ta.Ta) da = j {Vj^WtJITo) da 
Jo Jo 

= [ {VrSjiJlTa) da- [ {R{Ta,Ji)Ji,Ta) da 
Jo Jo 

= / Ta{Wjeji,Ta) da- f {R{Ta, 4) 4, Ta) da 

Jo Jo 
= [{"^J^JM]]- I {R{Ta,4)Ji,Ta) da 



{R{Ta,4)Ji,Ta) da 
using the fact that V jeJ^ — for a — 0, 1. So finally, 
dpt{x) = Y,(4{t,x),et{x)) dBl 



l>0 



+ ( Jo{u{t, gt{x)), u{t, gtix))),et{x) 



1 



e>o 

il.N I 



2pt(x) 

with J^-'^{t,x) the part of J^{t,x) normal to Tq 
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Removing the assumption that yt does not belong to the cutlocus of xt , it is well 
known (see [5] for a similar argument) that the formula becomes 



dpt{x) = ^ ( jo(i, x), et{x)) dB, 



+ {Jo{u{t, gt{x)),u{t, gtix))), etix)) ~ dLt{x) 



2pt{x) 



where ~Lt{x) is the local time of pt{x) when {gt{x), gt{x)) visits the cutlocus. Then 
letting 



Pt = p{gt,gt) = ( / Pt{x)dx 

IM 



1/2 



we get 



dpt^—y^i I Pt{x) (jo{t,x),et{x)) dx] dBf 
Pt ^ \Jm ^ I J 



+ — / Pt{x)[jii{u{gt{x)),u{gt{x))),et{x))dxdt~— [ pt{x)Lt{x)dx 
Pt Jm \ / Pt Jm 



+ 



/ V(4(t,a;),et(x)\ dx dt 
^Pt JMf^^\ I 

[ Pt{x){jo[t,x),et[x)^ dx\ dt. 

£^0 



For a vector w € Tg^^^^M, we let w'^ the part of w tangential to To{t, x). Letting 



dx dt 



(4'^(t,-),To(i,-) 



Im (4'^(^.2;),ro(t,a;)^ dx 



dx 



1/2 



(observe P( = / ||To(t, a;) ||^ c?x), we finally proved 
Jm 
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Proposition 6.1. The ltd differential of the distance pt between gt and gt is given 
by 



c^Pt = -7 X! ( / ^*(^) iP~3t{x),gt{x){(^I{9t{x))) - a-J{gt{x)) dx] dBf 
Pt \Jm / 

+ — / Pt{x){Pgt(a:),gtix){u^{gt{x)))) -u^{gt{x))) dxdt- — [ pt{x)dLt{x) dx 
Pt Jm Pt Jm 



pt Jm 



dxdt. 



In the case of manifolds with negative curvature we may observe a similar phe- 
nomena to the one of the torus with the Euclidean distance treated in Section 

4: as long as the norm stays sufficiently small to avoid the cut-locus of the 
manifold, the mean distance between the stochastic particles tends to increase 
exponentially fast. 



6.2. The rotation process. 

In the following we would like to study the rotation of two particles gt{x) and 
gt{x) when they are in a close distance one to another. Recall that we have noted 

Xt = gt(x), Ut = gt{x). We always assume that the distance from xt to yt is small: 
we are interested in the behaviour of e{t) as pt{x) goes to 0. We let 



(6.6) dmx{t)^ = a{xt)dBt - {a{xt)dBt, e{t))e{t) 



and 



(6.7) dmyit)'' = <T{yt)dBt - {a{yt)dBt, Pa^„yAt)) Px„yAt) 



where Pxt.-yait) denotes the parallel transport along 7a. 
^Prom Ito formula we have 



(6.8) ^To = pt{x)S>e{t) + dpt{x)e{t) + dpt{x)S>e{t) 
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and this yields 



^e(t) = -^^To h-dpt{x)e{t) - ]--^dpt{x)9e{t) 

Pt{x) pt{x) ^Pt{x) 

- ^ -jo{dmX{tf,dmy{tf) 



Ptix) 

{u{t, yt)) - u{t, xt), e{t)) e{t) 



Pt{x) 



2pt 

1 1 

2Mx 



-dpt{x)3ie{t) 



^ Md„,x{tf ,d^y{tf) + -^jo{u^{t,xt),u''{t,yt)) 



Pt{x) pt{x) 



2pt(a;) 



2p 



^ ( fT.i\\'^TXf-R{Ta,Ji)Ji,Ta) da] e{t) 



where we used the fact that dpt{x)^e{t) = 0, and where denotes the part of u 
which is normal to the geodesic 7o. Now as before 

V^.V^.To = VtoV^. Jo^ - i?(To, 4)4. 

Finally we get 
Lemma 6.2. 

^e{t) = -^jo{dmx{t)'',dmy{tf) + J—jo{u''{t,Xt),u''{t,yt)) 

Pt{x) pt{x) 



+ E Vto Vj. - i?(To, 4)4 dt 



t>0 



2p 



^ i tT.i\\'^Tjir-R{T.,Ji)4,Ta) da] eit). 



^Prom now on we assume that M = T the two dimensional torus. 

In this situation the curvature tensor vanishes and we have the formulas 

JaiVjW) = V + a{w — v), Ja{v,w) = W — V. 

We immediately get 

de{t) = 9e{t) (dmyit)'' - drnxit)"") + {{u''{t,yt) - u''{t,xt)) dt 

Pt[x) Pl(x) 
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where we used the fact that VtoVjcJ = 0, as a consequence of Vjf Jg = 0, 
Vj^^jf = 0, and R = 0. 

Let us speciahze again to the case where the vector fields are given by 



AkiO) = ik2,-ki)cosk.e, Bk{e) = (fc2,-fci)sinfc.0 
and the Brownian motion 

(6.9) dW{t) = J2 \kV^{Akdxk + Bkdyk) 

feez 

where Xk,yk are independent copies of real Brownian motions. As in section 2] we 
assume that J2k I^P-^fc < o° consider = A(|fc|) to be nonzero for a equal 

number of ki and ^2 components. Again we write 

(6.10) dgt = {odW{t)) + u{t, gt)dt, dgt = {odW{t)) + u{t, gt)dt 
with 

90 = .90 = ■0, </> 7^ ■0- 

Changing the notation to gt — gt {0) — xt, gt = gt {0) — yt, get 

de(t) = — \— Afe (cos k ■ gt — cos k ■ gt) k-^'^dxk 

H TTTT AfcVi' (sin k ■ gt - iwk ■ gt) k^'^ dyu 

+ ^Ai^''itr9t)-u^{t,gt)) dt 

- „ 2in\ (^(cosfc ■ gt-cosk- gtf + (sin A; ■ gt - sink ■ gt)^^ e(t) dt 

P*^ ' \k\^0 



1^ XkV^k^^^ (2sin^i:i^pM) dzk 



Pt\t>) 



J2xl.\k^^^\'sin'{^^\^)eit)dt 



^ ' |fe|#0 



where Zk is the Brownian motion defined by 



, . k-{gt+gt), , k-{gt+gt) 
dzk = — sm dxk + cos dyk- 
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Noting = \kf{nk ■ e(^))^ we obtain 

= ^ E \k\>^kMnk ■ e{t))e'{t) (2sm ^'\~ dzk 

(6.11) + ^ {{u^'itrgt) - u^it,gt)) dt 

Pt{0) 

^ J2 \k\'>^lHnk ■ e{t)f sin^ M^LzMe(i) dt 



where e'{t) is a unit vector in T orthonormal to e{t). Now for every i^T > 0, if 
7r 

Pt{d) < -^j^ then for all k such that |fc| < K, 

sin2 tlii^ 1 



Now using |fc| = \k^\ and {uk ■ e{t))^ + {n^.± ■ e{t))^ — 1, we get 

(6-12) E ■ Bin^ ^ ' %~ ''^ > ^ E ^'l^l'" 

^* ^ ' |fc|5:^0 a<\k\<K 

Observe that the term in the left is the second part of the drift in equation (|6.1ip 
as well as the derivative of the quadratic variation of e{t). This yields the following 
result. 

Proposition 6.3. Identifying TT with C, we have e{t) = e^'^^ where Xt is a real- 
valued semimartingale with quadratic variation 



.13) d[X, X], = ^Y1 \k\'>^Mnk ■ e(t)f sin^ ^J^^ dt 

' Ifcl^o 

and drift 



* 1 



(6.14) / —-—{u{s,gs)-u{s,gs),ie{s))ds. 
Jo Ps{0) 

We have for all K > 0, on ^Pf{9) < 

(6.15) d[X,X]t>^ J2 ^fcl^l'- 

Q<\k\<K 

If Xl\k\^ = +00, then as gt{6) gets closer and closer to gt{0), the rotation e{t) 
\k\^o 

becomes more and more irregular in the sense that the derivative of the quadratic 
variation of Xt tends to infinity. 
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